Efficiency and Its Bounds for a Quantum Einstein Engine at Maximum Power 
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We study a quantum thermal engine model for which the heat transfer law is determined by Ein- 
stein's theory of radiation. The working substance of the quantum engine is assumed to be a two- 
level quantum systems of which the constituent particles obey Maxwell-Boltzmann(M.B.), Fermi- 
Dirac(F.D.) or Bose-einstein(B.E.) distributions respectively at equilibrium. The thermal efficiency and 
its bounds at maximum power of these models are derived and discussed in the long and short ther- 
mal contact time limits. The similarity and difference between these models are discussed. We also 
compare the efficiency bounds of this quantum thermal engine to those of its classical counterpart. 
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INTRODUCTION 



The Carnot engine plays a crucial role in the founda- 
tions of thermodynamics. However, it can not be real- 
ized in practice since its output power is infinitesimally 
small due to its reversibility. Real thermal engines can 
not work as slowly as Carnot engines to preserve equi- 
librium and must lose energy during the working cycles 
due to various reasons. This makes the efficiency of real 
thermal engines below that of an ideal Carnot engine. 
To optimize the thermal engines in the real world, a lof 
of "realistic" models have been established and studied 
in the literaturejlrla]. One of the most practical problem 
associated with the optimization of real heat engines is 
its efficiency at maximum power. This problem was 
firstly studied by Curzon and Ahlborn in 1975ll|]. For 
Carnot-like heat engines, the authors assumed that the 
temperature differences between the heat reservoirs and 
working substance are finite and fixed, thus the two heat 
transferring processes are not reversible anymore, while 
the adiabatic expansion and compression processes are 
still reversible. Under these assumptions, they derived 
the well-known CA efficiency rycA = 1 — y/T c /Th, where 
Th and T c are the temperatures of the hot and cold heat 
reservoirs with which the working substance is in con- 
tact. Though the CA formula has a good agreement 
with measured efficiencies of some thermal plants, this 
model still has some intrinsic drawbacks. On one hand, 
it gives neither an exact nor constraint result for the ef- 
ficiency as pointed out by Ref.[4]. On the other hand, 
in real world situations the temperature differences be- 
tween the working medium and heat reservoirs are not 
constant and the heat transferring process could be gov- 
erned by some more general physical laws which can 
incorporate temperature changing during heat transfer- 



ring processes. 

In our previous work on classical engines|8], we have 
seen that the heat transfer law plays a crucial role on 
the efficiency at maximum power problem. For Carnot 
engines, the time periods for which the adiabatic expan- 
sion and compression processes last are usually negligi- 
bly short, while those of the two isothermal heat trans- 
ferring processes are infinitely long, therefore Carnot en- 
gines have zero output power. In real world situations, 
the isothermal heat transferring processes must last for a 
finite period of time and obey specific heat transferring 
laws. In Ref.|8], we studied a thermal engine model for 
which Newton's cooling law is obeyed during the heat 
transferring processes, and derived the upper and lower 
bounds for the efficiency at maximum power in the long 
and short contact time limits respectively. By consid- 
ering the heat transferring processes during which the 
temperature of the working medium is close to or far 
from isothermal, and adjusting the ratio between the 
heat capacities of the heating and cooling stages, the 
model can simulate different types of engines including 
but not limited to Carnot engines. 

The studies of classical thermal engines can be gener- 
alized to their quantum counterparts. Recently, differ- 
ent models of quantum thermal engines are extensively 
studied in the literaturep- lllll . The efficiency of a quan- 
tum thermal engine at maximum power has also been 
studied in Ref.[12], where the quantum thermal engine 
is based on the model discussed in Ref. 11311 , in which the 
quantum thermal engine is composed of particles con- 
fined in a one-dimensional (ID) infinite potential well, 
and the wall of the well can expand to perform work. 
The derived efficiency at maximum power is a univer- 
sal number. 

Following the same spirit, we try to generalize our 
previous work to the quantum world. For simplicity 
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and without losing generality, the working substance in 
our model of thermal engines is a two-level quantum 
system fl3ll , which can be chosen as the lowest two lev- 
els of a ID infinite quantum well or a ID harmonic oscil- 
lator. Unlike the discussion in Ref.Qjl], the temperature 
rather than average energy is used to describe the ther- 
mal equilibrium state as in Ref.[9]. What is important 
here is that the heat transferring process between the 
working substance and the heat reservoir is described 
by Einstein's theory of radiation. This can be thought as 
the quantum version of the model discussed in Ref.JH] in 
some way. For simplicity, we denote this kind of quan- 
tum thermal engine by "quantum Einstein engine". We 
are interested in the efficiency and its bounds at maxi- 
mum power. 

The organization of the paper is as follows. We first 
study the quantum thermal engine of which the con- 
stituent particles of the working substance obey the 
M.B.(Maxwell-Boltzmann) distribution. This will shed 
light on our successive studies on the other two mod- 
els. For this model, we derive the heat transferring law 
based on Einstein's theory of radiation, and give the for- 
mulas of the heat transfer and entropy production. As 
in Ref.|8j], we study the efficiency at maximum power 
and its bounds in the long and short contact time limits. 
We also study quantum engines for which F.D.(Fermi- 
Dirac) and B.E.(Bose-Einstein) distributions are applied. 



QUANTUM EINSTEIN ENGINE ASSOCIATED WITH 
M.B. DISTRIBUTION 



General results for heat transfer and entropy production 

The working substance of our quantum thermal en- 
gine is assumed to be a two-level quantum system with 
energy levels E\ (low) and E 2 (high). The energy differ- 
ence of the two levels is E 2 — E\ = hv, where h is the 
Planck's constant. The particle numbers at low and high 
energy levels are N% and N 2 respectively, and the fixed 
total particle number is given by Nq = Ni + N 2 . For 
simplicity, we first consider the case that the constituent 
particles of the quantum system satisfy M.B. distribu- 
tion at equilibrium states. Assume the initial tempera- 
ture of the system is T\, then the initial particle distribu- 
tions are 



Ni =N 



1 



1 + cxp (—fiihv) 



,N 2 = N 



1 



1 + cxp (Pihv) 



, (1) 



where $\ = 1/fceTi and ks is the Boltzmann constant. 
In our model of quantum thermal engine, the heat reser- 
voir can be thought of as a black-body source with tem- 
perature T 2 . When the working substance or the two- 
level quantum system is "in contact with" a black-body 
source, the heat is transferred by the photon emission 



and absorbtion. This heat transferring process is de- 
scribed by Einstein's theory of radiation[14] 



dN 2 



dt dt 

= BN Q u v - 2BN 2 u u - AN 2 = -aN 2 + b, (2) 

where a = 2Bu u + A, b = BN f) u u . A and B are the 
famous Einstein's coefficients, and u u is the spectral en- 
ergy density of the black-body source. The solution of 
Eq.© is given by 

N 2 (t) = --[-- N 2 (0)] exp (-at) (3) 
a a 

Introducing the distribution function /(/3) = 1/(1 + 
exp (f3hv)) for level E 2 , then we have 

iVa(O) = Nof(px), N 2 {t -> 00) = - = N f(/3 2 ). (4) 

a 

At time t, we have N 2 (t) = f((3(t))N where 

f((3(t)) = /(&) - [/(&) - f(fii)] exp (-at). (5) 

Similarly, for level E\, we have Ni(t) = (1 - f(J3(t)))N . 
We are interested in the situation with j3hv <C 1, which 
is true if the size of the quantum well is not too small or 
the spring constant of the harmonic oscillator is not too 
large at high temperature. To the leading order of (3hv 
we get 

1 1 94 

/(/3) ~ - 7/^, a = Acoth(p 2 hv/2)n—. (6) 

Therefore, Eq.01 leads to the changing of inverse tem- 
perature 

2At , 



P(t)=h~(p2~Pi)eM 



hv$ 2 



(7) 



It is interesting to notice that the time constant 1 /a de- 
pends on the temperature of the heat reservoir. If tem- 
perature is higher, less time is needed for the work- 
ing substance to reach equilibrium with the heat reser- 
voir, this phenomenon is counterintuitive as shown by 
Fi gur e IH Note that in general Einstein's coefficient A oc 
v 2 [15], then the inverse time constant a is in fact propor- 
tional to v. Therefore the higher the energy difference 
hv, the longer the time is needed for the system to reach 
the equilibrium state. 

In terms of the distribution function /(/?), the quan- 
tum entropy of the two-level system can be expressed 
as 

S = -k[f(p) ln/(/3) + (1 - f(p)) In (1 - /(/?))]. (8) 

Plug in the expression of /(/3), and expand Eq.© in se- 
ries of fihv 



S = k\n2 



-k(fihv) 2 + 0(/3hv) 



(9) 
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FIG. 1: (Color online) Temperature of the working substance 
during the heat transferring process as a function of time. The 
blue solid (red dashed) line indicates the temperature chang- 
ing of the working substance for the situation that the tem- 
perature (T 2 ) of the heat reservoir is 400K (1000K). For both 
situations, the initial temperature of the working substance is 
300K. The unit is chosen such that A = 1 and hv/k — 1. 



For a heat transferring process between the initial state 
with inverse temperature /3o and the final state with in- 
verse temperature (3(t) at time t, the heat transfer and 
entropy production are given by 

AQ(t) = (N 2 (t) - N 20 )hv = N hv(f(p(t)) ~ f(M) 

' 1,2,2 



= -N h 2 v 2 (f3 -f3(t))+O((t3hp) 3 ), 



AS(t) = -fc(H 2 (/3 2 -/3 2 W) + O((/3H 4 ), 



(10) 



where N20 and N 2 (t) are the particle numbers of initial 
and final state respectively. In the derivation of AQ(t), 
we have used the fact dN 2 — —dN%. Generically, the 
energy levels of the quantum working substance dur- 
ing the heating and cooling stages are different since the 
engine must expand to perform work between the two 
stages. We assume that the frequencies associated with 
the heating and cooling stages are Vh and v c respectively 
and the initial and final inverse temperatures are fih§, 
Ph{t) and /3 c0 , /3 c (t) respectively. Therefore, to the lead- 
ing order of fihv the heat transfer and entropy produc- 
tion during the heating stage are given by 



AQ h (t) = -N a h 2 v 2 h (l3 m ~ &(*)), 



Similarly during the cooling stage we have 



AQ c {t) = -N h 2 v 2 ((3 c0 - (3 c (t)), 



-kh 2 v 2 Mz-m\ 



(12) 



We assume the time duration that the heating and cool- 
ing stages last are and r c respectively. When the quan- 
tum engine finishes a full thermodynamical cycle, the 
working medium returns back to its initial state and we 
have ASh{Th) + AS c (t c ) — 0. The power output and the 
efficiency of the thermal engine are given by 



P 
V 



AQh(r h ) + AQ c (t c ) 

Th + T c 

1 | AQ c (r c ) 
&Qh(Th)' 



(13) 



Here we adopt the convention that Q > 0(< 0) means 
absorbing(releasing) heat. As mentioned before, the 
working substance of the quantum engine under con- 
sidering is composed by particles at the two lowest lev- 
els of a ID infinite quantum well or ID harmonic os- 
cillator. The cycles of the quantum engine are as fol- 
lows. l)The system absorbs heat Q/, from the hot reser- 
voir during the time period r/j . 2)The system expands to 
perform work. In case of ID infinite quantum well, the 
width of the well expands from Lq to L\. In case of ID 
harmonic oscillator, the"width" of the harmonic poten- 
tial expands while its spring constant "shrinks" from kg 
to fci in its parameter space. 3)The system releases heat 
Q c at the cold reservoir during the time period r c . 4)The 
system returns to its original size and temperature and 
the working medium returns to its original state. 



Efficiency and its bounds in the long contact time limit 

When the contact time is long enough, the working 
substance can exchange heat sufficiently with the reser- 
voirs. In this limit, r^. c / '(f3hvh, c A) is large, and we as- 
sume the final inverse temperatures for the two stages 
are j3h = Ph{Th) arid j3 c — /3 c (r c ). Hence the heat trans- 
fers during the two stages are 

AQ ft = \N h 2 v 2 x, AQ C - -lN h 2 p 2 c y, (14) 

where x = ftho — fih and y = f3 c — /3 c q. To the leading 
order of fihv, the constraint ASh + AS C = gives 



v 2 h {2p h +x)x-v 2 c (2p c -y)y = 0, 
of which the solution is given by 



AS h (i) 



(ii) 



0* + —(20 c -y)y-p h , 



(15) 



(16) 
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where 7 = Vh/v c . Substitute Eq,l fT6|l to the expression of 
the output power 



4 T h +T c 



(17) 



and let dP/dy = 0, the only meaningful solution for y is 
found to be 



y= TT7 2 (18) 

Plug y into Eq. fTSl we obtain the thermal efficiency at 
maximum power 



Vn 



= 1-^ 



(19) 



1 



7 2 (1 -y c )-l + y/(l + 7 2 )[1 + 7 2 (1 - Vc) 2 
(2-^ c )7 2 



where rj c = 1 — ^ is the Carnot efficiency. If we let 
7 approach and 00 respectively we obtain the upper 
and lower bounds of r\ m as 



Vc . , ??c 



(20) 



Interestingly these bounds agree exactly with those 
given in Ref.J3] for classical thermal engines in the 
long contact time limit. However, the situation is a 
little different for our model of quantum thermal en- 
gine. The quantum engine, either quantum-well type or 
harmonic-oscillator type, must expand to perform work 
after absorbing heat at hot reservoir. Hence its size must 
increase afterwards. For a quantum-well type engine, 
its width L will increase, while for a harmonic-oscillator 
type engine, its spring constant k will decrease. Since 
the frequency is anti-proportional to L 2 or proportional 
to \fk, then 7 = t>h/v c must be larger than 1, or the lower 
limit of 7 is 1 rather than 0. Now we have tighter bounds 
for r\ m 



2 - ^4 - 4t? c + 2tj2 



< Vm < 



Vc 



2 - 77 c 2 — 77 c 

The efficiency r\ m can be expanded in series of 77 c as 



1 



'In 



7% 



7 



1 1 o(ug). 



(22) 



2' ,u ' 4(1 + 7 2 ) 

The coefficient of the second-order term lies between 1 /8 
and 1/4. 

Efficiency and its bounds in the short contact time limit 

In the short time limit such that r <C I /a, the inverse 
temperature can be approximated to the second order of 
ar as 

(3(t) a / 3 1 + ( / 3 2 _ / 3 1 ) ar + ( /3l _ / 3 2 )I a V 2 . (23) 



Implementing this approximation, the entropy produc- 
tions during the two stages are given by 

AS h = isl[-2/3 h0 (l3 h - /3 h0 )a h T h 

-(PhO-PhKWhO-Malr 2 }, 
AS C = i/ c 2 [-2/3 c o(/3 c -/3 c o)a c r c 



-(/3 c0 -/3 c )(2/3 c0 -/3 c )a 2 



(24) 



Using the same convention as in the last subsection, the 
constraint ASh + AS C = gives an equation for x and y. 
Since a^, c r^, c is an infinitesimal quantity, we match both 
sides of the equation order by order of ah, c Th,c- To the 
first and second order, we get 

7 2 (/^ + x)x-^-r h = (f3 c - y)y-^r c , 

PhVh PcVc 

1 2 x(p h + 2x)a 2 h rl = (fi c - 2y)ya 2 c r 2 , (25) 



from which one can deduce 

Ph + 2x 



x(/3 h + xf 



= 7 



v{pc - yf 



(26) 



The key step here is to simplify the above equation and 
get a relatively simple relation between x and y as in 
Ref.[8], hence we can avoid messing up the physics 
by the mathematical complexity. We assume that the 
temperature difference is small relative to the tempera- 
ture of the heat reservoir at each heat transferring stage. 
Thus, x is small relative to /3h and y small to f3 c . Expand- 
ing both sides of Eq. d26b to the third order of x and y, we 
can derive a very simple relation between x and y 



x B, 



(27) 



The heat transfers during the two heating and cooling 
stages are given by 

AQ h = ^N h 2 v 2 h xa h T h , AQ C = -^N h 2 v 2 ya c T c ,(28) 



(21) 

from which the output power is given by 



P = -^N h 



1 , 2 v 2 h xa h T h ~ v 2 c ya c T c 



Th 



(29) 



Taking into account that the spontaneous emission coef- 
ficient A satisfies A <x v 2 , except a constant the output 
power is evaluated as 



P (X 



(& - i 2 xB h /B c -B h - x)x 
B c - i 2 xB h /B c + 7 (B h + x)B 2 JB 2 ' 



(30) 



Solve the equation dPj dx = 0, and plug the only mean- 
ingful solution for x into the expression of the efficiency 
at maximum power 



B h + x 



Bc-l 2 xB h /8 c - 



(31) 
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In series of T) e , r\ m can be expanded as 



^7m — 



Vc , 7(27 2 +7 + l) 



-»£ + 0fa c 3 ). (32) 



2 ' 8(l + 7 + 7 2 +7 3 )' 
Taking 7 = 0, oo, the bounds of J] m are derived as 

-77 <Vm< 7, • (33) 

2 2 - rj c 

Interestingly, we obtain the same rough bounds of r\ m as 
in the long time contact limit. Again noting that 7 can 
not be smaller than 1, the finer bounds of f] m are found 
to be 



1 - yjl-ric < T] m < - 

2 - ilc 



(34) 



The coefficient of the second-order term of r\ m also lies 
between 1/8 and 1/4. 



The quantum entropy of the working substance is ap- 
proximated by 



ln2-— p 2 h 2 v 2 . 
32^ 



(38) 



One can see that the expressions of f ¥ (f3), a F and 5 F dif- 
fer from their maxwellian counterparts only in the coef- 
ficients of the leading order of j3hu. If we go on carrying 
the calculations as before, it is not difficult to find that 
the results are the same as those associated with M.B. 
distribution. In other word, the quantum Einstein en- 
gine with fermionic working substance has no signifi- 
cant difference from that with the maxwellian working 
substance if j3E\ 2 <C 1. 



QUANTUM EINSTEIN ENGINE ASSOCIATED WITH 
THE F.D. DISTRIBUTION 



QUANTUM EINSTEIN ENGINE ASSOCIATED WITH 
THE B.E. DISTRIBUTION 



In this section, we consider the case that the con- 
stituent particles of the working substance obey the F.D 
distribution. In fact we will see that this situation will re- 
duce to that associated with the M.B. distribution even 
when /3£i j2 < 1 (If /3£i, 2 > 1, both ED. and B.E. distri- 
butions reduce to M.B. distribution). Here we still adopt 
the same convention for the parameters as in the last sec- 
tion. For each level of the quantum system, the particle 
numbers of the initial state are given by 



iVi(O) 



N Q 



e Pl E 2+l 



N 2 (0) 



N 



(35) 



Similar to what we have done in the last section, the 

fFfa\ - 1 I ft i e^+i- 



is m- 



distribution function f — J./V- 1 ~<~ e ^ E i+i 
troduced. Thus the initial number distributions are 
iVi(O) = JV (1 - fWi)) and N 2 (0) = N o f(0i). When 
the quantum system is in contact with the hot (cold) 
reservoir, the heating (cooling) process is described by 
Einstein's theory of radiation Eq.©. Solve this equa- 
tion, we get Nx(t) = N (l - f(j3(t))) and N 2 (t) = 
N f ¥ ((3(t)) in which f v (f3(t)) is also expressed by Eq.|0 
with (3(0) — /?i and /3(oo) = (3 2 . Hence the heat ab- 
sorbed by the working substance at time t is AQ(t) = 
Nohv(f(J3{t)) - f((3(0))). To the leading order of 0hv 
we have 



f((3) « i - Iphv, 



^2 -Ei 



4A 



(36) 
(37) 



General results for heat transfer and entropy production 

Now we consider the the quantum Einstein engine of 
which the working substance obeys the B.E. distribu- 
tion. The discussions follow quite similar steps as the 
previous section. The initial particle distributions of the 
two-level quantum systems are 



iVi(O) = 



N, 







1 



e' j i E i -1 
e' J i E 2-l 



JV 2 (0) = 



N, 







e/ 9lE 2 _l 



(39) 



wrw; 



The heating (cooling) process is again governed by Ein- 
stein's theory of radiation. Introducing the distribution 
function / B (/3) = 1/(1 + ff^f ), then at time t the par- 
ticle distributions become N t (t) = N (l - f B (/3(t))) and 
N 2 (t) = N f(f3(t)) with 13(0) = fa and /3(oo) = /3 2 . 
The heat transfer is given by AQ(t) = Nohis(f B (/3(t)) - 
/ B (/3(0))) . Similarly, to the leading order of f3hv we have 



1 hvf3E x E 2 



Ei + E2 2 {E\ 
A(E 1 + E 2 ) 



(40) 



The quantum entropy production is further approxi- 
mated by 



6 



Ei 



■In- 



Ei 



+ - 



In- 



£2 



.E1 + E2 Ei + £2 -Ei + -E2 -Ei + E2 



E1E2I11 fafti/ EiE 2 hv(3(E$ - El) + 2{Ef - AE1E2 + El) In fi) 



2(Ei + E 2 ) 2 



24(Ei + E 2 ) 3 



J41) 



In what follows, we will ignore the superscript "B" for 
simplicity. For the heating stage, we use i^o and Eh to 
denote the initial (low) and final (high) energy levels, 
Pho and f3h (t) to denote the initial and final inverse tem- 
peratures. Define the frequency vk by Eh — Eho = hvh, 
then the heat transfer and entropy production during 
the heating and cooling stages are given by 

AQ h (t) = \N h 2 v 2 h X h (p h0 - p h {i)), 

AS h (t) = Y hl (p h0 -p h (t))+Y h2 (p 2 h0 -p 2 h (t)),(42) 

AQ c (t) = ^N h 2 is 2 X c ((3 c0 -p c {t)), (43) 

AS c (t) = Y cl (f3 c0 ~f3 c (t))+Y c2 (f3 2 c0 -f3 2 (t)) : (44) 
where 



Yh: 



Y c2 



(E h0 + E h y 



E h E h0 ln^-hv h 
2(E h + E h0 ) 2 ' 



E h0 E h hu h (3(El - E* ) + 2(E 2 - 4E h0 E h + E 2 ) In 

24(E h0 + E h )3 ' 
E c0 E c _ E c E c0 ln^hu c 

' (■ ] 



(E c0 + E c ) 2 ' " c± 2(E C + E c0 ) 2 
E c0 E c hu c (3(E 2 - E 2 ) + 2(E 2 - 4E c0 E h + E 2 C ) In §^ 
2A(E c0 + E c ) 3 



It is maximized when dP/dy = 0. Let •y 2 = jf^^r, then 
we have ^ — ^7 =0. Using Eq.|46|, the only meaningful 
solution for y is 



_ P'dil + 7 2 4 ) - V(P' 2 if + g_g +~ 2 ) uq . 

The efficiency of the thermal engine at maxim power is 



-1 Xcv 2 c y 1 

X h Vh x 

Plugging Eqs.([46t and J49] | we have 



(50) 



j 2 7l(l - rQ - 1 + J-^/(l + 7 ?(1 - jgjaj (j + 7 4) 

1 7 2 2 7 2 2 +7i 2 (l-^) 

(51) 

where 77^ = 1 — can be thought of as the corrected 
Carnot efficiency. If we chose suitable Eh, c and £*o,cO 
such that 71 = 72 = 7', Eq.lTSTb "recovers" the result (fl9)l 
except that r/ c is replaced by ^ and 7 by 7'. By inspecting 
the numerators and denominators of 71 and 72, one can 
find that they are of the same order of Eh c and Eho, c o 
respectively. Hence it is reasonable to assume that they 
are of the same order when approaching or 00. There- 
fore we get a rough estimation of the upper and lower 
bounds 



Efficiency and its bounds in the long contact time limit 

From now on, the discussions are simply parallel to 
what we have done in the last section. We briefly outline 
our results here. In this limit, Th. c —> 00. As previously 
did, we assume Phij h ) = Ph, Pc{t c ) = Pc, Pho - Ph = % 
and (3 c o — /3 C — —y. From ASh + AS C = we have 

Yhix + Y h2 (2p h + x)x - Y cl y - Y c2 {2p c - y)y = 0. (45) 

Solve this equation, we have 



x = Jp' h 2 + ^m-y)y~PL 



nl a 1 Yhl , Y c \ 2 Yh 2 

^^^ 2 ^^ty- 2 ^ = y7 2 

The power of this heat engine is given by 

1 .Ji^i - X c vly 



P : 



-Noh 2 ' 



Th 



(46) 



(47) 



(48) 



2 - Vm - 2 - rf c 



(52) 



Following the same reasoning as before, the tighter 
bounds of r] m are given by 



2 - y/4 - 4^ + 2^ 
2-V'c 



< Vm < 



(53) 



Despite the similarity between the results associated 
with the B.E. distribution and those associated with the 
M.B. and FD. distributions, there is also qualitative dif- 
ference between them. Obviously the latter only depend 
on the difference of energy levels, or Vh,c- However, the 
former has an explicit dependence on the choice of ini- 
tial (low) energy level Eho, c o- 



Efficiency and its bounds in the short contact time limit 

We consider the limit that r <C 1/a and expand the 
inverse temperature to the second order of or, then the 
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entropy productions during the two stages are given by 

AS h = -2Y h2 (/3' h + x)xa h T h -Y h2 (p' h + 2x)xali%, 
AS C = 2Y c2 {B' c -y) V a c T c + Y c2 {fi' c ~2y)ya 2 c T 2 . (54) 

While the heat transfers are 



1 



AQfc = -N h v h X h xa h T h , 



AQ C = --^N h v c X c ya c T c 



(55) 



As before we match both sides of ASh + AS C = order 
by order of ar, then the first and second orders of ar 
give 



7i (P'h + x)xa h T h = (f3' c - y)ya c T c , 



7 2 (f3 h + 2x)xa 2 T 2 = {p' c -2y)ya, 
from which we deduce 



2 T 2 



x(f3' h + x) 2 ll y(P' c ~y) 



P' c -2y 



(56) 



(57) 



Using the same argument as in the last section, one can 
find 



x 11 B' 



(58) 



The spontaneous emission coefficient A satisfies A cx v 2 , 
then the output power is evaluated as: 



{I3' c -llxl3' h /p' c -\{p' h +x))x 

TTW> (59) 



where 73 = ^rffi^^r ■ Th e power is maximized when 

dP/dx = 0, which leads to the solution of the efficiency 
at maximum power is 



fl'v*4./r-v2<v„ -vi / / 3 ^h(/ 3 c+^h7 2 2 )(l+73)(/3^;+^7 2 2 73) 
Phil + Pcl2li - 71 U p,2 + pg^ 

~?ra,^ R >^ / /9^1(1+73)(^7 1 2 +/3^ 2 2 7 3 ) 1 



(60) 



If Eh c and E^q.co are chosen suitably such that 7! = 
72 = 73 = 7'/ the efficiency 77™ can also be expanded in 
series of r]' c as 



(l+7' + 7 /2 +7 /3 ) 



■rg + 0(rg). (61) 



Interestingly this also "recovers" the result (|32l l expect 
that ?7 C is replaced by t]' c and 7 by 7'. Similarly when 7' 
approaches 0, 00 or 1, 00, the rough and fine bounds of 
rjm are found to be 



Vc ^ , ^ 

2 ~ Vm ~ 2~V C 



1 - V 1 - V'c < Vm < 



(62) 



CONCLUSIONS 



In conclusion, we presented an analysis of the quan- 
tum thermal engine of which the heat transferring law 
is derived from Einstein's theory of radiation. Again we 
notice that heat transferring laws play a crucial role on 
the problems about the efficiency of quantum thermal 
engines at maximum power. The thermal efficiency and 
its bounds at maximum power for quantum Einstein en- 
gines are studied in the long and short time limits. To 
some extent, this can be thought as the quantum coun- 
terpart of the classical thermal engine studied in Ref.[8], 
which can simulate some well-known classical thermal 
engines. For j3hv <C 1, We find that the quantum 
Einstein engine with fermionic working substance has 
no difference from that with maxwellian working sub- 
stance, while the one with bosonic working substance 
has a qualitative difference. 
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